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PREFACE 


This  report  was  prepared  by  the  Environmental  Laboratory  (EL)  of  the 
US  Army  Engineer  Waterways  Experiment  Station  (WES) ,  as  part  of  the  Environ¬ 
mental  Impact  Research  Program  (EIRP) ,  Work  Unit  No.  31730,  "Environmental 
Impacts  of  Modifying  Estuarine  Circulation  and  Transport  Processes."  The  EIRP 
is  sponsored  by  the  Office,  Chief  of  Engineers  (OCE) ,  US  Army,  Washington,  DC. 
The  OCE  Technical  Monitors  for  EIRP  are  Dr.  John  Bushman  and  Dr.  David  Buelow. 
Mr.  Dave  Mathis  is  the  Water  Resources  Support  Center  Technical  Monitor. 

This  report  describes  the  use  of  a  box-type  model  for  transport/water 
quality  calculations  and  the  linkage  of  this  model  to  multidimensional  hydro- 
dynamic  models.  Performance  of  the  coarser  grid/time-step  box  model  is  com¬ 
pared  with  directly  linked  transport  codes. 

The  study  was  conducted  and  the  report  was  prepared  by  Ms.  Sandra  L. 

Bird  and  Mr.  Ross  Hall  of  the  Water  Quality  Modeling  Group  (WQMG) ,  Ecosystem 
Research  and  Simulation  Division  (ERSD) ,  EL,  under  the  direct  supervision  of 
Mr.  Mark  Dortch,  Chief,  WQMG.  General  supervision  was  provided  by 
Mr.  Donald  L.  Robey,  Chief,  ERSD,  and  Dr.  John  Harrison,  Chief,  EL. 

Dr.  Roger  T.  Saucier  was  Program  Manager  of  EIRP.  The  report  was  edited  by 
Ms.  Lee  T.  Byrne  of  the  Information  Products  Division,  Information  Technology 
Laboratory. 

Commander  and  Director  of  WES  during  the  preparation  and  publication  of 
this  report  was  COL  Dwayne  G.  Lee,  CE.  Dr.  Robert  W.  Whalin  was  the  Technical 
Director. 

This  report  should  be  cited  as  follows: 

Bird,  Sandra  L.,  and  Hall,  Ross.  1988.  "Coupling  Hydrodynamics  to  a 

Multiple-Box  Water  Quality  Model,"  Technical  Report  EL-88-7,  US  Army 

Engineer  Waterways  Experiment  Station,  Vicksburg,  Miss. 
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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used 
(metric)  units  as  follows: 

_ Multiply _ 

cubic  feet 
degrees  (angle) 
feet 
inches 

pounds  (mass) 
square  feet 


in  this  report  can  be  converted 

gy 

0.02831685 

0.01745329 

0.3048 

2.54 

0.4535924 

0.09290304 


to  SI 

To  Obtain 
cubic  metres 
radians 
metres 
centimetres 
kilograms 
square  metres 
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COUPLING  HYDRODYNAMICS  TO  A  MULTIPLE-BOX  WATER  QUALITY  MODEL 


PART  I :  INTRODUCTION 

1.  Assessment  and  management  of  water  quality  are  aided  by  the  use  of 
numerical  models.  Water  body  hydrodynamics  interact  with  biological,  chemi¬ 
cal,  and  other  physical  processes  to  affect  water  quality  variables.  Many 
simplifying  assumptions  are  typically  made  in  the  equations  governing  these 
processes  in  order  to  develop  economical  numerical  models.  These  assumptions 
are  made  in  relationship  to  a  particular  situation  (i.e.,  an  assumption  that 
might  be  acceptable  for  analyzing  a  problem  in  one  water  body  might  provide 
fallacious  results  in  another). 

2.  One  of  the  most  routinely  used  assumptions  suppresses  variation  of 
the  variables  within  a  cross  section,  and  the  model  equations  are  written  in  a 
one-dimensional,  longitudinal  form.  Typically,  the  cross-sectionally  averaged 
approach  is  used  for  riverine  applications.  For  a  few  well-mixed  homogeneous 
estuaries,  this  approach  may  also  be  appropriate.  Reservoir  models  sometimes 
assume  no  variation  in  the  horizontal  dimension  and  solve  the  one-dimensional, 
vertical  form  of  the  equation.  One-dimensional  models  such  as  CE-QUAL-R1 ,  a 
vertical  reservoir  model  (Environmental  Lab  (EL)  1986)  and  CE-QUAL-RIV1 ,  a 
longitudinal  riverine  model  (EL,  in  preparation)  solve  the  transport/water 
quality  equations  quickly  and  efficiently.  The  one-dimensional  assumption, 
however,  limits  the  problems  that  can  be  adequately  addressed  with  these  types 
of  models.  Only  some  reservoir  problems  and  a  very  limited  number  of  problems 
in  estuaries  and  coastal  embayments  can  be  adequately  addressed  using  this 
one-dimensional  approach. 

3.  In  recent  years,  many  two-dimensional  and,  even  more  recently, 
three-dimensional  hydrodynamic  models  have  been  developed  and  applied  to  res¬ 
ervoirs,  estuaries,  and  coastal  embayments.  No  single  model  can  appropriately 
describe  currents  and  mixing  in  all  of  these  water  bodies.  Highly  stratified 
estuaries  require  consideration  of  vertical  variation  of  velocity  and  water 
quality  constituents;  wide  estuaries  require  consideration  of  lateral  varia¬ 
tions;  and  large  estuaries  (which  may  be  both  wide  and  stratified)  can  require 
resolution  in  all  three  spatial  dimensions.  Reservoirs  may  be  deep  and  strat¬ 
ified  or  broad  and  shallow.  Because  of  this  variety  in  water  bodies,  several 
different  two-  and  three-dimensional  hydrodynamic  models  have  been  developed 


U 


and  used  in  estuarine  and  reservoir  applications  at  the  US  Army  Engineer 
Waterways  Experiment  Station  (WES)  by  the  Hydraulics  Laboratory  (HL)  and  the 
Coastal  Engineering  Research  Center  (CERC). 

4.  The  necessity  of  evaluating  environmental  impacts  of  US  Army  Corps 
of  Engineer  (USACE)  activities  on  a  variety  of  water  bodies  was  the  impetus 
for  the  development  of  multidimensional  water  quality  modeling  capabilities  by 
the  EL.  Three  considerations  guided  the  selection  and  development  of  a  multi¬ 
dimensional  water  quality  modeling  approach: 

a.  Long-term  multidimensional  water  quality  modeling  can  become 
cumbersome  and  computationally  very  time  consuming  when  water 
quality  algorithms  are  directly  linked  to  hydrodynamic  models. 

b.  The  EL  should  be  able  to  perform  water  quality  modeling  studies 
in  conjunction  with  hydrodynamic  studies  performed  by  both  the 
HL  and  CERC. 

c.  Water  quality  kinetics  rarely  require  the  spatial  and  temporal 
resolution  required  for  accurate  hydrodynamic  calculations.  To 
meet  these  requirements,  a  multiple-box  (also  known  as  a  mixed 
segment,  cells  in  series,  or  integrated  compartment)  model  was 
chosen  as  the  transport  framework  for  a  versatile,  computa¬ 
tionally  efficient,  water  quality  model.  This  type  of  model  can 
be  overlaid  on  the  same  grid  as,  or  a  coarser  grid  than,  the 
hydrodynamic  model,  and  it  can  use  a  larger  time-step.  Hydro- 
dynamic  model  output  can  be  averaged  over  time  and  space  to 
drive  the  water  quality  model. 

5.  This  report  first  describes  the  formulation  and  limitations  of  a 
multiple-box  model.  The  Water  Quality  Analysis  Simulation  Program  (WASP), 
developed  under  the  auspices  of  the  US  Environmental  Protection  Agency 
(USEPA) ,  was  adapted  for  the  purposes  of  this  study;  and  the  linkage  of  WASP 
to  two  different  hydrodynamic  models  is  described.  Transport  applications  of 
the  multiple-box  model  were  made  to  three  different  water  bodies  with  differ¬ 
ent  physical  characteristics:  the  Savannah  River  Estuary,  DeGray  Reservoir, 
and  the  Mississippi  Sound.  Results  are  analyzed  and  discussed. 
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PART  II:  MULTIPLE-BOX  MODEL  FORMULATION 


General  Formulation 


6.  WASP,  the  USEPA  multiple-box  model,  was  adapted  for  use  in  this 
study.  Concentrations  in  this  model  are  determined  by  simple  mass  balance 
around  a  series  of  completely  mixed  reactors,  from  the  following  equation 
(Ditoro,  Fitzpatrick,  and  Thomann,  1983): 


vi  3T  '  SQJ,icJ  +  J  -  ci>  !  “i  !  Vi 


(1) 


(2) 


(3) 


(4)  (5) 


(1) 


where 


^  =»  segment  volume ,  L 

i  =  segment  index 

j  =  index  of  adjoining  segment 


J » i 


=  net  advective  flow  from  segment  j  to  segment  i  ,  L  /T 


=  concentration  in  segment  j  ,  M/L 
E^  ^  *  dispersion  coefficient  for  the  i,j 


cross-sectional  area  of  the  i,j 


interface,  L  /T 

2 

interface,  L 


=  mixing  length  between  segment  i  and  j 


i.J 
Li.j 

=  segment  concentration,  M/LJ 

W  «  point  or  distributed  sources  and  sinks  of  the  constituent,  M/T 

1  3 

=  kinetic  degradation  or  transformation  rate,  M/L  T 

This  ordinary  differential  equation  is  solved  in  the  WASP  model  using  Euler's 

Method. 


7.  Figure  1  is  a  schematic  illustrating  the  processes  described  by 
Equation  1.  "A  "  an<*  "^out"  (term  2  in  Equation  1)  represent  the  flux  of 
material  into  and  out  of  the  cell  by  the  net  advection  of  the  velocity  field. 
"D"  (term  3  in  Equation  1)  represents  the  concentration  flux  from  some  equal 
exchange  flow  between  the  two  cells  accounting  for  dispersive  transport 
resulting  from  velocity  and  concentration  fluctuations  across  the  dimensions 
averaged.  "B"  (term  4  in  Equation  1)  is  the  flux  of  material  to/from  the 
segment  boundary.  "K"  (term  5  in  Equation  1)  represents  the  change  in 
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concentration  arising  from  reactions  occurring  within  the  segment  including 
both  degradation  and  transformation  reactions. 

8.  The  multiple-box  model  is  the  result  of  a  volume  average  for  each 
segment  of  the  three-dimensional  advective  diffusion  equation;  i.e.,  it 
becomes  a  zeroth-dimensional  type  equation.  Individual  box  segments  can  then 
be  arranged  in  any  arbitrary  manner  forming  a  one-,  two-,  or  three-dimensional 
network.  Although  Figure  1  is  an  example  of  a  one-dimensional  alignment  of 
segments,  the  processes  shown  can  occur  between  two-  and  three-dimensional 
segment  arrangements  as  well.  This  geometric  flexibility  allows  the  segments 
to  intermesh  with  any  hydrodynamic  model  grid  as  long  as  box  volumes  and  flows 
between  the  boxes  can  be  calculated. 

Dispersive  Properties 

9.  Dispersion  presents  the  most  difficulty  in  the  application  of  the 
multiple-box  model.  Two  major  problems  arise  regarding  dispersion  in  the 
application  of  the  multiple-box  model.  First,  calculation  of  an  appropriate 
dispersion  coefficient  for  use  in  the  multiple-box  model  is  difficult. 

Second,  the  multiple-box  model  may  be  numerically  overdif fusive .  The  numer¬ 
ical  diffusion  introduced  by  the  solution  technique  may  be  greater  than  the 
physical  dispersion  of  the  system  unless  relatively  small  segments  are  used, 

10.  Calculating  dispersion  for  any  transport  model  is  difficult.  In  a 
dimensionally  averaged  model,  the  primary  contribution  to  dispersion  arises 
from  nonuniformity  of  concentration  and  velocity  in  the  dimension  of  aver¬ 
aging.  Although  the  choice  of  dispersion  coefficients  is  very  difficult  for 
one-  and  two-dimensional  estuarine  transport,  some  systematic  guidelines  are 
available  in  the  literature  (Fischer  1976  and  Fischer  et  al.  1979).  For  the 
multiple-box  model,  no  systematic  guidelines  are  available  in  the  literature 
for  estimating  these  parameters. 

11.  However,  one  potential  procedure  for  adjustment  of  the  dispersion 
coefficient  is  based  on  duplicating  dye  transport  simulated  with  the  hydrody¬ 
namic  model.  Typically,  if  the  multiple-box  model  is  used  in  conjunction  with 
a  multidimensional  hydrodynamic  model,  the  hydrodynamic  model  will  include  a 
transport  algorithm  for  calculating  salinity  in  an  estuary  or  temperature  in  a 
reservoir.  Hydrodynamic/transport  models  are  generally  calibrated  against 
field  measurements  of  salinity  and/or  temperature  distribution  or  a  dye  study 
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Figure  1.  Schematic  representation  of  box  model  processes 


in  the  water  body.  Injection  and  transport  of  a  conservative  tracer  in  the 
hydrodynamic  model,  compared  with  the  transport  of  an  identical  injection  in 
the  multiple-box  model  over  a  short  period  of  time  (i.e.,  a  period  of  time 
typical  for  the  hydrodynamic  model  runs)  ,  can  be  used  as  a  guide  for  the 
adjustment  of  the  dispersion  coefficient  in  the  multiple-box  model.  Addi¬ 
tionally,  this  procedure  would  be  effective  in  assessing  the  errors  introduced 
in  the  use  of  tidally  averaged  values  for  advective  transport. 

12.  Excessive  numerical  diffusion  is  often  a  critical  limitation  in  the 
multiple-box  concept.  Shanahan  and  Harleman  (1984)  evaluated  the  diffusive 
properties  of  a  multiple-box  model  using  a  one-dimensional  arrangement  of 
boxes,  compared  with  a  one-dimensional  advective  dispersion  model  for  steady 
uniform  flow 


+  Qi£  =  Ai_  ^ 
at  ^ax  ax  ax 


(2) 


where 

•> 

A  =  cross-sectional  flow  area,  M*- 

3 

Q  =  system  through  flow,  M  /T 

x  =  the  coordinate  in  the  direction  of  flow,  L 

2 

D  =  one-dimensional  dispersion  coefficient,  M  /T 
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Their  analysis  is  based  on  conceptual  reactor  models  used  in  sewage  treatment 
plant  design.  A  peclet  number  (Pe)  is  defined  as 


Pe 


QX 

AD 


(3) 


where  X  is  the  total  length  of  the  system  in 
analysis  indicates  that,  to  keep  the  box  model 
diffusive,  there  must  be  at  least  n  segments 
being  a  large  number) 


the  direction  of  flow.  Their 
from  being  inherently  over- 
where  n  is  defined  as  (for 


n 


n 


(4) 


They  point  out  that  for  a  given  spatial  step  size.  Ax  =  X/n  ,  Equation  U  is 

equivalent  to  numerical  diffusion  in  an  upwind  spatial  finite  difference 

approximation  of  Equation  2,  i.e.,  D  =  l!  x/2  ,  where  D  is  the  numerical 

n  n 

diffusion  coefficient  introduced  by  the  solution  scheme  and  IJ  is  the  cross- 
sectionallv  averaged  velocity  il./T)  . 

13.  However,  according  to  Roache  (1982),  numerical  diffusion  for  upwind 
differencing  in  a  one-dimensional  system  is  of  the  form 

T!  Av 

D  =  i-H(!  -  u)  (5) 

n  ^ 


where 


EAt 

X 

At  =  time 


a  is  referred  to  as  the  Courant  number.  As  the  Courant  number  approaches  one, 
numerical  diffusion  approaches  zero.  The  condition  where  a  =  1  is  the 
stability  limit  for  the  upwind  differencing  scheme.  In  a  nne-dimefisionai 
system  when  a  =  1  ,  then  the  value  of  At  represents  the  time  that  it  takes 
for  a  particle  to  travel  the  length  of  a  cell  (Ax).  This  stability  criterion 
applied  to  the  box  model  takes  the  form  At/V^  -  1  and  can  be  interpreted 


as  the  total  flow  into  or  out  of  a  segment  during  a  time-step  must  not  exceed 
the  volume  of  the  segment.  This  interpretation  lends  itself  to  extrapolation 
to  multidimensional  problems.  Restating  Equation  5  in  terms  of  multiple-box 
model  parameters,  numerical  diffusion  at  each  segment  interface  in  the 

multiple-box  model  can  be  described  for  equal  length  segments  by 


D 

n 


Q.  .1.  . 

.-\ll  1 

2Ai,  j 


(6) 


One-Dimensional  Example 

14.  The  transport  properties  of  the  multiple-box  model  are  illustrated  by 
considering  steady  uniform  flow  in  a  rectangular  channel  with  the  following 
characteristics: 

H  =  1.0  ft* 

W  =  20.0  ft 
L  =  2,000  ft 
U  =  0.2  fps 

Mn  =  0.017 
D  =  14 . 7  f c‘/sec 
A  =  20.0  ft^ 

Q  =  4.0  c  f  s 

Pe  =  27.2,  peclet  number 

where  H  is  the  channel  depth,  W  is  the  channel  width,  L  is  the  total 
channel  length,  U  is  the  average  longitudinal  velocity,  and  Mn  is 
Manning's  roughness  coefficient.  The  dispersion  coefficient,  D  ,  was  cal¬ 
culated  using  the  method  of  Fischer  et  al.  (1979).  For  a  very  small  time-step 
a  <<  1  ,  the  channel  must  be  divided  into  at  least  14  segments  (n  >  Pe/2)  to 
avoid  excessive  numerical  diffusion  in  the  model. 


*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  presented  on  page  3. 
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15.  The  effect  of  increasing  the  time-step  on  numerical  diffusion  in 
the  box  model  is  illustrated  in  Figure  2.  For  an  initial  concentration  of 
10  mg/ l  in  the  most  upstream  box,  concentration  versus  time  profiles  in  box 

13  calculated  with  four  time-steps  ranging  from  0.001  to  0.0082  days  are  shown 
in  this  figure.  No  physical  dispersion  was  input  for  these  test  cases. 
Numerical  diffusion  decreases  as  the  time-step  increases,  dropping  dramati¬ 
cally  as  the  time-step  limit  (At  =  0.008275  days)  is  approached.  This  behav¬ 
ior  is  consistent  with  Equation  6. 

16.  The  box  model  results  were  compared  with  the  analytical  solution  of 
the  one-dimensional  advective  diffusion  equation  for  an  instantaneous  point 
source  of  material  injected  into  steady  flow  in  a  uniform  channel.  The 
solution  is  given  by  Crank  (1984): 


C(x, t) 


2 

-x 

4Dt 


(7) 


where  x  is  the  downstream  distance  from  the  point  of  injection,  t  is  the 
time  since  the  injection,  M  is  the  mass  injected,  and  D  is  the  one¬ 
dimensional  dispersion  coefficient.  Figure  3  shows  the  results  for  an  injec¬ 
tion  of  M  =  809  g  (the  amount  of  mass  equal  to  10  mg/ i  injected  into 

segment  1  of  the  14-segment  box  model  discretization)  injected  at  x  *  0  , 

2 

t  -  0  ,  and  a  value  of  14.7  ft  /sec  for  the  dispersion  coefficient  D  com¬ 
pared  with  a  14-segment  box  model  simulation  using  a  0.001-day  time-step.  The 

2 

numerical  diffusion  for  the  box  model  simulation  was  calculated  as  12.6  ft  / 

2 

sec.  A  value  of  2.1  ft  /sec  was  input  for  the  dispersion  coefficient  in  the 

2 

model  to  give  a  total  model  dispersion  of  14.7  ft  /sec.  Peak  concentrations 
are  slightly  lower  in  the  upstream  segments  of  the  multiple-box  model  since 
the  initial  mass  injection  must  be  spread  over  the  entire  box,  rather  than 
being  a  true  point  source  injection.  This  difference  is  reduced  at  the  end  of 
the  channel,  and  box  model  transport  approaches  the  analytical  solution. 


Model  Modifications 


17.  WASP  was  developed  for  lake  applications  using  constant  volume 
boxes  and  steady  flows.  The  constant  volume  and  steady  flow  assumptions  are 
not  acceptable  for  intratidal  estuarine  applications;  i.e.,  the  time-step  is 
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CONCENTRATION,  MG; 


SEGMENT  NUMBER  3 


SEGMENT  NUMBER  6 


Figure  3.  Comparison  of  box  model  solution  to  analytical 
solution  for  flow  in  a  one-dimensional  channel 


less  than  a  tidal  cycle.  For  variable  volume  applications,  the  volume  on  the 
left  hand  side  of  Equation  1  must  be  written  within  the  time  differential, 
i.e.,  D(V^C^)/dt  .  The  value  for  the  quantity  was  found  using  a  Euler 

solution  scheme  and  solving  for  the  concentration  at  the  new  step.  The 
modified  formulation  conserves  mass  for  both  constant  and  variable  volume 
applications.  For  variable  volume/unsteady  flow  applications,  the  model  must 
be  altered  to  read  values  for  volumes  and  flows  at  every  computational  time 
rather  than  to  read  them  once  during  the  initial  data  input. 
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PART  III:  HYDRODYNAMIC  INTERFACING 


General 


18.  In  this  study,  the  multiple-box  model  is  interfaced  with  hydrody¬ 
namic  output  generated  by  two  different  models.  The  first  is  CE-QUAL-W2,  a 
two-dimensional,  laterally  averaged  hydrodynamic  model  developed  for  the  USACE 
(EL  and  HL  1986).  This  model  was  originally  developed  for  two-dimensional 
reservoir  modeling  and  extended  for  use  in  deep,  narrow,  stratified  estuaries 
with  the  addition  of  estuarine  boundary  conditions  (Edinger  and  Buchak  1981). 
Applications  of  the  Savannah  River  Estuary  (Hall  1987)  and  DeGray  Lake  (Martin 
1987)  are  used  herein  as  case  studies  for  utilization  of  CE-QUAL-W2  output  as 
the  hydrodynamic  driver  of  a  multiple-box  model. 

19.  The  second  model  used  to  generate  hydrodynamic  output  for  the 
multiple-box  model  in  this  study  is  WIFM-SAL  (WES  Implicit  Flooding  Model  with 
constituent  transport)  (Schmalz  1985b)  ,  a  vertically  averaged  model  employing 
an  exponentially  stretched  grid.  WIFM-SAL  was  developed  for  the  analysis  of 
shallow  estuaries  and  embayments  that  could  be  assumed  to  be  vertically  well 
mixed.  An  application  of  WIFM-SAL  to  the  Mississippi  Sound  and  adjacent  areas 
(Schmalz  1985a)  was  used  as  the  test  case  for  the  interface  with  the  multiple- 
box  model. 


Interface  with  CE-QUAL-W2 


Savannah  River  application 

20.  Figure  4  shows  the  computational  grid  for  the  main  channel  and  tide 
gate  branches  for  the  application  of  CE-QUAL-W2  to  the  Savannah  River  Estuary 
with  the  multiple-box  model  segments  overlaid  on  it.  A  relatively  coarse  grid 
overlay  of  18  box  model  segments  was  made  on  a  total  of  377  active  computa¬ 
tional  cells  in  CE-QUAL-W2.  In  the  upper  end  of  the  estuary  where  the  reach 
is  predominantly  riverine  and  unstratified,  a  single  vertical  layer  of  boxes 
was  overlaid  on  the  CE-QUAL-W2  grid  and  expanded  to  a  double  layer  in  the  par¬ 
tially  stratified  downstream  sections.  Thus,  fine-scale  vertical  resolution 
was  not  maintained  in  this  box  model  overlay. 

21.  Variables  in  CE-QUAL-W2  are  defined  as  shown  in  Figure  5.  Water 
surface  elevation  (Z) ,  cell  width  (B) ,  and  constituent  concentrations  are 
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GEOMETRIC  SCHEMATIZATION 
SAVANNAH  RIVER  ESTUARY  APPLICATION 
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Figure  4.  Savannah  River  Estuary  grid  for  CE-QUAL-W2  with 
box  model  overlay  shown  by  heavier  lines 
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Figure  5.  Variable  definition  sketch  for  CE-QUAL-W2 


defined  at  the  center  of  each  cell.  Layer  thickness  is  a  constant  (H)  for  all 

except  the  top  layer,  which  is  variable  (H-Z) .  Velocities  (u,w)  and  diffusion 

coefficients  (D  ,  D  )  are  defined  at  the  cell  faces, 
x  z 

22.  Box  model  segment  volumes  were  calculated  at  each  CE-QUAL-W2  time- 
step  by  summing  the  CE-QUAL-W2  cell  volumes  within  each  box  segment.  The 
CE-QUAL-W2  cell  volume  fV(I,K> ]  was  calculated  by 

V  ( I ,  K.)  »=  B(I,K)  x  H(K)  x  Ax  K  >  KT 


V(1 ,K)  =  B(I,K)  x  [H(K)  -  Z ( I ) ]  x  Ax 


K  =  KT 


(8) 


The  flow  across  each  box  segment  face  was  calculated  by  summing  corresponding 
CE-QUAL-W2  cell  flows  at  each  box  face.  K  is  the  layer  number,  and  KT  is 
the  surface  layer.  CE-QUAL-W2  flows  are  given  by 


Qh(I  -  l.K)  =  U(I  -  1  ,K)  x  h(K)  x  B(I-’K)  1  K  ^  KT,  I  >  1 

Qh(I  -  2  ,KT)  =  U(I  -  1  ,KT)  x  .Z2(I)  lH^T.)  -  Z2(I  -  1)] 

x  -+  »»•*>!  K  -  KT,  !  >  1 


(I ,K  -  1)  =  W(I,K  -  l)  x  Ax  x 


[B(I,K  -  1)  +  B(1 ,K) ] 
2 


K  >  1 


(9) 


where  Qh  and  represent  flows  in  the  horizontal  and  vertical  directions. 

The  flow  across  the  upstream  face  at  1=1  into  the  first  multiple-box 
segment  is  simply  the  upstream  boundary  flow  specified  in  CE-QUAL-W2 .  Box 
model  flows  are  calculated  by  summing  the  flows  across  each  of  the  CE-QIJAL-W2 
cell  faces  that  align  with  the  box  model  face. 

23.  The  time-averaged  flow  (Q^)  was  calculated  as  the  arithmetic  aver¬ 
age  of  the  flows,  as  follows: 


N  Q 

=  y  — 

^ ,  N 


n=l 


(10) 


where  N  is  the  number  of  time-steps  averaged.  However,  in  the  time-averaged 
data  set,  the  volume  at  the  beginning  of  each  averaging  interval  was  used 
since  the  average  net  flow  into  a  segment  over  the  averaging  interval  added  to 
the  volume  at  the  beginning  of  the  interval  equalled  the  volume  at  the  begin- 
ning  of  the  next  interval.  In  this  way,  continuity  was  assured. 

DeGray  Lake  application 

24.  Whereas  the  Savannah  River  is  a  strongly  advective  system  with  a 
residence  time  on  the  order  of  days,  DeGray  Lake  typically  has  a  residence 
time  of  several  months  and  exhibits  very  strong  thermal  stratification  during 
the  summer.  Relatively  fine  vertical  resolution  is  required  for  accurate 
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water  quality  modeling.  The  DeGray  application  of  CE-QUAL-W2  used  a  horizon¬ 
tal  segment  length  of  993.6  m  and  a  layer  thickness  of  2.0  m.  The  first  box 
model  test  used  a  direct  grid  overlay  resulting  in  428  segments  and  828  flows. 
The  correspondence  between  CE-QUAL-W2  cells  and  box  model  cells  in  the  first 
box  model  test  is  shown  in  Table  1. 

25.  A  second  box  model  test  consisted  of  a  2  by  2  overlay  on  the 
CE-QUAL-W2  grid,  i.e.,  four  CE-QUAL-W2  cells  per  box  model  segment  except  for 
some  of  the  bottom  segments.  This  box  model  application  consisted  of  112  seg¬ 
ments  and  210  flows.  Correspondence  between  CE-QUAL-W2  cells  and  box  model 
cells  in  the  second  box  model  test  is  shown  in  Table  2.  Daily  averaged  values 
for  flow  and  volume  input  for  both  of  these  box  model  overlays  were  calculated 
using  the  same  general  approach  as  previously  described  for  the  Savannah  River 
application. 

26.  The  CE-QUAL-W2  time-step  for  the  DeGray  application  was  1,500  sec 
or  0.02+  days.  The  box  model  time-step  was  0.1  day.  The  1,500-sec  time-step 
size  selected  for  CE-QUAL-W2  reflects  an  internal  gravity  wave  restriction.  A 
subset  of  CE-QUAL-W2  hydrodynamic  output  calculated  during  the  summer  months 
was  used  for  model  comparisons,  since  interest  centered  primarily  on  erosion 
of  the  thermocline  resulting  from  numerical  diffusion  in  the  box  model. 
Experimentation  revealed  that  0.1-day  time-steps  were  computationally  stable 
for  the  box  model  during  the  time  interval  documented  in  this  report;  however, 
during  autumnal  overturn,  the  time-step  size  in  the  box  model  was  limited  to 
0.02  day  because  of  the  Courant  number  restriction.  Daily  averaged  values  for 
the  volumes  and  flows  were  used  repeatedly  for  the  time-steps  in  a  particular 
simulation  day. 


Interface  with  WIFM-SAL 


27.  Figure  6  shows  the  transformed  coordinate  grid  with  the  locations 
of  specific  variables  used  in  WIFM-SAL  calculations.  The  velocities  (U,V)  are 
defined  at  the  cell  faces  while  depth  (h)  ,  water  surface  elevation  (n) ,  and 
constituent  concentration  are  defined  at  the  cell  center.  The  variables 
and  a9  are  the  spatia1  coordinates  in  transformed  space.  Figure  7  shows  the 
computational  grid  in  real  space  coordinates  for  the  WIFM  application  to  the 
Mississippi  Sound  and  surrounding  areas.  The  WIFM  grid  is  59  by  115,  i.e., 
nearly  7,000  cells.  More  boxes  were  needed  to  provide  an  adequate  overlay  on 
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WASP-01  Overlay  on  CE-OUAL-W2  Grid* 
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WASP-01  Overlay  on  CE-QUAI.-W2  Grid* 
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Figure  6,  Variable  definition  sketch  for  WIFM 

this  system  than  were  required  in  the  Savannah  River  application.  The  over¬ 
lay,  selected  to  provide  adequate  resolution  within  the  constraints  of  reason¬ 
able  storage  requirement  and  computation  time,  was  a  regular  2-  by  3-overlay 
on  the  WIFM  grid  (bold  outlines  on  Figure  7),  i.e,,  6:1  WIFM  cells  per  box 
model  segment  resulting  in  nearly  1,000  box  model  segments.  The  use  of  a 
regular  overlay  made  it  possible  to  automate  generation  of  the  interface  file 
that  defined  box  model  segments  and  flow  faces  in  terms  of  the  WIFM  cells. 

This  information  is  required  for  the  generation  of  the  box  model  input  from 
hydrodynamic  input.  For  a  large  number  of  box  model  cells,  manual  generation 
of  this  interface  is  tedious  and  time  consuming. 

28.  An  approach  slightly  different  from  that  described  for  CE-QUAL-W2 
was  used  in  calculating  volumes  and  flows  for  box  model  input  from  WIFM 
hydrodynamics.  The  solution  scheme  for  the  continuity  equation  employed  in 
WIFM  was  used  as  a  basis  for  calculation  of  multiple-box  volumes  and  flows. 

An  approach  analogous  to  that  of  Schmalz  (1985b)  in  the  development  of  a 
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Sound  computational  grid  for  WTFM-SAL  with  1:6  box  model 
over  lav  shown  with  heavy  lines 


three-time-level  explicit  transport  scheme  for  use  in  WIFM-SAL  was  employed. 
The  finite  difference  forms  of  the  continuity  equation  used  in  the  alternating 
difference  type  solution  of  the  two-dimensional  hydrodynamics  were  combined  to 
yield  one  three-time-level  finite  difference  expression. 
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where 


nn".n  ~  water  surface  elevation  at  time  level  kil  in  cell  (n,m) 
it  =  time-step  length 

(p.)  =  stretching  coefficient  in  a  direction  at  cell  index  n 

i  m  L 

Au ,  =  ^  direction  space  increment 

uk_i  ,  =  velocity  component  in  a,  direction  lor  ceil  (n,m)  at  tin 
n’Q+;  level  kil  1 

( U  , ) n  =  stretching  coefficient  in  direction  at  cell  index  m 


v  _  =  velocity  component  in  direction  for  cell  (n,m)  at  time 

n  ’ ‘n  level  kii 


d  -  water  depth  in  cell  (n,m)  at  time  level  k 
n ,  m 

29.  Rearranging  Equation  11  and  adding  the  time  invariant  depth  -h 
to  the  water  surface  elevation  yields 


where  term  1  is  a  finite  difference  expression  for  volume  of  cell  (n,m)  at  the 
k+1  time  level,  which  is  set  equal  to  the  volume  of  the  cell  at  the  k-1 
time-step  (term  2) ,  plus  the  net  volume  change  resulting  from  the 
approximations  for  flow  over  two  time  Intervals  across  the  four  faces  of  the 
cell  (terms  3,  4,  5,  and  6).  Continuity  can  be  guaranteed  in  the  generation 
of  box  model  parameters  if  calculation  of  volume  and  flow  are  based  on 
Equation  12. 

30.  Box  segment  volume  is  calculated  as  the  sum  of  the  volume  of  the 

k-1 

cells  overlaid  where  individual  WIFM  cell  volume,  V’ol  ,  is  calculated  (as 

n  ,m 

suggested  by  Equation  12)  by 


k-1 

n  -  h 
n )  m  n  ,m 


(  V, )  Aa  (  y„)  Aa 
1  m  !  l  n  l 


Likewise,  box  model  flows  are  calculated  by  summing  the  flows  across  each  of 
the  WIFM  cell  faces  aligned  with  the  box  model  face.  The  flow  into  each  WIFM 
cell  is  found  for  rx^  and  directions  by  dividing  terms  3  and  5  respec¬ 

tively  by  2  At  (terms  4  and  6  represent  flow  out  of  the  WIFM  cell).  Since 
Equation  12  is  a  three-time-level  finite  difference  expression,  flows  and 
volumes  are  calculated  at  alternate  WIFM  time-steps  since  terms  3  and  5  give 
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the  flow  from  the  k-1  to  k+1  time  levels  and  the  time  interval  used  for 
multiple-box  input  is  2At  where  At  is  the  WIFM-SAL  time-step.  Additional 
time-averaging  can  then  be  performed  as  described  for  the  interface  with 
CE-QUAL-W2. 
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PART  IV:  TRANSPORT  APPLICATIONS 


General 


31.  To  identify  the  potential  and  limitations  for  simulating  multi¬ 
dimensional  transport  using  the  relatively  coarse  grid,  long  time-step, 
multiple-box  model,  movement  of  a  conservative  constituent  in  the  box  model 
was  compared  with  movement  of  the  same  constituent  in  the  finer  scale, 
directly  linked,  transport  models.  Results  for  the  box  model  and  CE-QUAL-W2 
simulations  were  compared  using  both  Savannah  River  Estuary  and  DeGray  Lake 
applications;  comparison  with  WIFM-SAL  used  the  Mississippi  Sound  application. 

Comparison  with  CE-QUAL-W2 


Savannah  River  application 

32.  A  25-mg/J.  instantaneous  dye  injection  was  made  in  segment  3  of  the 

box  model  and  in  the  equivalent  area  of  the  CE-QUAL-W2  grid  for  the  Savannah 
River  Estuary,  as  shown  in  Figure  4.  CE-QUAL-W2  hydrodynamics  were  averaged 
as  described  in  Part  III  and  used  to  drive  the  box  model.  Simulations  were 
performed  for  3.5  days.  CE-QUAL-W2  simulations  required  a  2-min  time-step, 
whereas  the  box  model  was  run  at  a  series  of  time-steps  up  to  3.5  hr.  In 
order  to  evaluate  the  ability  of  the  box  model  overlay  to  replicate  CE-QUAL-W2 
transport,  a  set  of  three  graphical  displays  was  made  for  several  of  the  box 
model  segments.  The  first  graph  was  a  time-history  of  the  volume  weighted 
average  of  the  concentrations  in  the  CE-QUAL-W2  cells  contained  within  a  given 
box  model  segment  (represented  in  the  graph  by  a  solid  line)  and  the  range  of 
the  concentrations  (shown  with  the  vertical  bar  (I))  found  in  these  cells 
(Figures  8a-14a).  The  second  graph  in  the  set  was  a  time-history  of  concen¬ 
trations  for  the  3.5-hr  box  model  simulations  (o _ o)  compared  with  the 

volume  weighted  CE-QUAL-W2  results  (A _ A)  (Figures  8b-14b) .  The  third 

graph  in  the  set  was  concentration  time-histories  in  the  segment  for  the  box 
model  simulation  using  different  time-step  sizes  (Figures  8c-14c)  ,  repre¬ 
senting  0.5-hr  (A _ A),  2-hr  (+ _ +) ,  and  3.5-hr  (o _ o)  time-step 

results,  respectively. 

33.  The  CE-QUAL-W2  cells  that  overlay  segment  4  immediately  downstream 
of  the  injection  location  show  the  greatest  variation  in  concentration  for  any 
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of  the  segments  (Figure  8a) .  The  coarse  box  model  overlay  did  not  provide  the 
resolution  of  the  finer  CE-QUAL-W2  grid.  Excessive  information  was  lost  where 
the  variation  of  concentration  within  the  segment  was  large  compared  with  the 
average  concentration  in  the  segment.  However,  the  average  concentration  in 
the  CE-QUAL-W2  cells  and  the  box  model  simulation  agree  (Figure  8b)  except  for 
the  concentration  peak  value  in  this  segment.  The  peak  concentration  in  the 
box  model  simulation  was  sensitive  to  the  time-step  choice  (Figure  8c).  For 
example,  the  0.5-hr  time-step  for  the  box  model  simulation  substantially 
decreased  the  peak  concentration  compared  with  the  2.0-  and  3.5-hr  time-steps. 
At  segment  5,  a  large  variation  in  the  concentrations  in  the  overlaid 
CE-QUAL-W2  cells  diminished  substantially  downstream  (Figures  9a-12a)  at  seg¬ 
ments  6,  7,  and  9.  The  box  model  replicated  the  oscillation  in  concentrations 
caused  by  the  flow  reversals.  The  phase  as  well  as  the  magnitude  of  the 
oscillations  was  generally  matched. 


Figure  8.  Savannah  River  application:  concentration  time-history 
for  box  model  segment  U  and  overlaid  CE-QUAL-W2  cells.  (a)  Vol¬ 
ume  weighted  average  of  CE-QUAL-W2  cells  overlaid  by  the  box  model 

segment  ( _ )  and  the  range  of  concentration  in  these  cells 

(|),  (b)  comparison  of  box  model  (o _ o)  and  CE-QUAL-W2 

(A _ A)  results,  (c)  comparison  of  box  model  results  using  a 

0.5-hr  (A _ A),  2-hr  (+ _ +)  ,  and  3.5-hr  (o _ o)  time-step 
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Figure  9.  Savannah  River  application:  concentration  time-history 
for  box  model  segment  5  and  overlaid  CE-QUAL-W2  cells. 

(a)  Volume  weighted  average  of  CE-QUAL-W2  cells  overlaid  by  the 

box  model  segment  ( _ )  and  the  range  of  concentration  in 

these  cells  (|),  (b)  comparison  of  box  model  (o _ o)  and 

CE-QUAL-W2  (A _ A)  results,  (c)  comparison  of  box  model 

results  using  a  0.5-hr  (A _ A),  2-hr  (+_ _ +) ,  and  3.5-hr 

(o _ o)  time-step 
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Figure  10.  Savannah  River  application:  concentration  time-history 
for  box  model  segment  6  and  overlaid  CE-QUAL-W2  cells.  (a)  Volume 
weighted  average  of  CE-QUAL-W2  cells  overlaid  by  the  box  model 

segment  ( _ )  and  the  range  of  concentration  in  these  cells 

(|),  (b)  comparison  of  box  model  (o _ o)  and  CE-QUAL-W2 

(A _ A)  results,  (c)  comparison  of  box  model  results  using  a 

0.5-hr  (A _ A),  2-hr  (+ _ +)  ,  and  3.5-hr  (o _ o)  time-step 
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Figure  11.  Savannah  River  application:  concentration  time-history 
for  box  model  segment  7  and  overlaid  CE-QUAL-W2  cells.  (a)  Volume 
weighted  average  of  CE-QUAL-W2  cells  overlaid  by  the  box  model 

segment  ( _ )  and  the  range  of  concentration  in  these  cells 

((),  (b)  comparison  of  box  model  (o _ o)  and  CE-QUAL-W2 

(A _ A)  results,  (c)  comparison  of  box  model  results  using  a 

0.5-hr  (A _ A),  2-hr  (+_ _ +)  ,  and  3.5-hr  (o _ o)  time-step 
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Figure  12.  Savannah  River  application:  concentration  time-history 
for  box  model  segment  9  and  overlaid  CE-QUAL-W2  cells.  (a)  Volume 
weighted  average  of  CE-QUAL-W2  cells  overlaid  by  the  box  model 

segment  ( _ )  and  the  range  of  concentration  in  these  cells 

(|),  (b)  comparison  of  box  model  (o _ o)  and  CE-QUAL-W2 

(A _ A)  results,  (c)  comparison  of  box  model  results  using  a 

0.5-hr  (A _ A),  2-hr  (+ _ +)  ,  and  3.5-hr  (o _ o)  time-step 
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TIM*  DAYS 

Figure  13.  Savannah  River  application:  concentration  time-history 
for  box  model  segment  11  and  overlaid  CE-QUAL-W2  cells. 

(a)  Volume  weighted  average  of  CE-QUAL-W2  cells  overlaid  by  the 

box  model  segment  ( _ )  and  the  range  of  concentration  in 

these  cells  ( | ) ,  (b)  comparison  of  box  model  (o _ o)  and 

CE-QUAL-W2  (A _ A)  results,  (c)  comparison  of  box  model 

results  using  a  0.5-hr  (A _ A),  2-hr  (+ _ +) ,  and  3.5-hr 

(o _ o)  time-step 


10  r 


SEGMENT  NO  13 


(a) 


0  0  5  in  15  2  0  2  5  in  35 


SEGMENT  NO  13 


SEGMENT  NO  13 


(c) 


Figure  1A.  Savannah  River  application:  concentration  time-history 
for  box  model  segment  13  and  overlaid  CE-QUAL-W2  cells.  (a)  Vol¬ 
ume  weighted  average  of  CE-QUAL-W2  cells  overlaid  by  the  box  model 

segment  ( _ )  and  the  range  of  concentration  in  these  cells 

(|),  (b)  comparison  of  box  model  (o _ o)  and  CE-QUAL-W2 

(A _ A)  results,  (c)  comparison  of  box  model  results  using  a 

0.5-hr  (A _ A),  2-hr  (+ _ +) ,  and  3.5-hr  (o _ o)  time-step 
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34.  However,  in  segment  11  (Figure  13b)  the  box  model  does  not  ade¬ 
quately  simulate  the  average  concentration  in  the  CE-Ql'AL-W'2  cells  that  it 
overlays.  The  box  model  underestimates  the  peak  concentration  in  this  segment 
even  using  the  3.5-hr  time-step.  The  source  of  this  error  becomes  clear  when 
the  large  range  of  concentrations  in  the  CE-QUAL-W2  cells  overlaid  by  seg¬ 
ment  4  (Figure  8a),  the  segment  upstream  from  segment  11,  Is  considered. 
Circulation  in  the  estuary  moved  the  highest  concentration  material  into  the 
bottom  cells  of  this  segment,  but  the  box  model  transported  the  average  value 
from  segment  4  into  segment  11.  Civen  the  particular  box  model  overlay  and 
injection  condition,  this  anomaly  was  the  result  of  lowered  spatial  resolution 
of  the  box  model  overlay.  At  segment  13  (Figure  14)  farther  downstream  in  the 
bottom  layer  of  boxes,  the  box  model  simulation  mimicked  the  average 
CE-QFAL-W2  results.  The  box  overlay  can  impact  accuracy  of  results  and  should 
be  carefully  considered  in  terms  of  the  problems  addressed.  Tf  a  sharp  front 
is  not  simulated,  results  will  not  be  so  severe. 

DeCray  Lake  application 

35.  The  primary  objective  of  the  DeGrav  hake  application  was  to  compare 
the  vertical  spreading  of  material  in  the  WASP  grids  with  the  CE-QUAL-W2 
simulation.  A  uniform  injection  of  dye  was  made  in  the  top  layers  of  both  the 
box  model  and  CE-QUAL-Vl'  grid.  For  the  1:1  WASP  (referred  to  as  WASPOI) 
overlay,  a  100-mg/f  iniection  was  made  into  the  surface  layer  of  cells  in  both 
models.  For  the  2:1  WASP  overlay  (referreu  tc  as  WASP02) ,  the  mass  injected 
into  the  surface  layer  of  the  CE-0LAL-W2  cells  was  distributed  through  the 
respective  WASP  surface  layer  segments.  The  initial  WASP02  segment  tracer 
concentrations  are  listed  in  Table  3. 

5b.  Daily  averaged  values  of  flows  were  calculated  from  the  CF-QhAL-W? 
output  for  WASP  input;  a  time-step  of  0.!  dav  was  used  in  the  WASP  simulation 
(i.e.,  each,  set  of  averaged  hydrodynamics  was  used  for  10  time-steps).  For 
the  WASPOI  simulation,  0.1  day  approached  the  time-step  limit;  for  the  WASP02 
simulation,  a  0.5-day  time-step  could  be  used  without  instabilities.  However, 
all  results  presented  here  used  a  0.1-day  time-step  in  the  simulation.  Jf 
simulations  were  carried  through  fall  overturn,  shor.ter  time-steps  had  to  be 
used,  as  mentioned  in  Part  III. 

37.  Concentration  versus  elevation  plots  at  six  longitudinal  segments 
alter  30,  60,  and  90  days  of  simulation  are  shown  respectivelv  in  Fig¬ 
ures  15-17.  In  these  figures,  TFST01  refers  to  the  concentration  in  the 


32 


tfh 


u 


co’Otf: 


C'j  ’LJ  I 

Us) 


LJ 
2 
oO 
O  LJ 


O 


nr; 


ci 
tr: 
o  *— 


OiJ'09  00*01’ 

1  1 U  A  3  3  3 


L;  a  x 


Oo 

OLJ 


LJ 

C 

cc 


O0’O0f 


l-Hrl-I-M-M-H - f -  * 

oo*oo(  oj*f»  cj'oy 
IVS)  NO  1 1 U AJ 33 


c>l 

3 

I 


oO 

OCJ 


0C  *0c ' 


U:  Q  X 


- - -  -■  y  - - i 

\iM  -•] rrry  oo  ‘u 

( ^ ;  NO  1 1 U A  333 


i-  _  -- 
»  1  \ ,f 
s;  £  3 

W  £0  X 


LJ 

>  vi 


O  u 
,.o 


— - i  - -  J 

ocroei  oo  *  oo  i  ■.:■  c-o*o9' 

thM  NO  I  luAj"-3 


X... 


S  * 
Sax 


oo  ■  ae i 


••Uj 


ooooi  un 'Oh 
(Wl  NOI1UA313 


z 

O  LJ 

Sec 

UJ 

L_> 

LX 

a: 

o^ 

o 


oo’of 


z&ti 


□  X 


l - - 

oo'oei  oo’ooi  o j ■ on 

IW)  N01ldA313 


UJ 

LJ 


00'0^ 


34 


FJuure  13,  DeCrny  Lake  L’ASJ'O;  test  refill  tf  compared  with  CF-QVA1 
results  after  30-day  simuJatfon 
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Figure  16.  DeGray  Lake  WASP01  test  results  compared  with 
CE-QUAL-W2  results  after  60-day  simulation 
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Figure  I7.  DeGrav  Fake  WASP01  test  results  compared  with 
('.H-QUAL-W2  results  after  90-day  simul  ion 


overlaid  CE-Q11AI.-W2  cells.  Results  were  very  similar  for  the  WASP01  and 
TEST01  simulations.  U'ASPOl  was  slightly  less  vertically  diflusive  than 
TEST01.  Numerical  diffusion  may  be  slightly  lower  for  WASP01  because  of  the 
larger  time-step.  While  the  CE-QUAL-W2  simulation  required  over  2  hr  of  simu¬ 
lation  time,  the  WASP01  required  less  than  4  min. 

38.  Similar  plots  were  prepared  for  the  WASP02  simulation  (Fig¬ 
ures  18-20).  In  these  plots,  the  TEST02  values  refer  to  the  volume  weighted 
average  of  the  concentrations  in  the  CE-QUAL-W2  cells  overlaid  by  the  WASP 
segments.  After  30  days  of  simulation  (Figure  18),  concentration  differences 
in  the  upper  layers  were  significant.  The  greatest  deviation  between  the  two 
simulations  occurred  in  the  upper  layers.  A  large  portion  of  this  deviation 
was  probably  due  to  the  fact  that  the  initial  box  concentrations  represented  a 
distribution  of  surface  CE-QUAL-W2  cell  tracer  mass  over  two  layers.  By 

60  days  (Figure  19),  the  tracer  was  generally  mixed  through  the  epilimnion  in 
both  cases,  and  the  only  remaining  significant  deviation  was  in  the  metalim- 
nion  at  the  segment  near  the  dam.  This  trend  continued,  and  at  90  days,  the 
deviation,  even  in  the  downstream  segment,  was  further  reduced.  The  increased 
concentration  observed  for  WASP02  in  segment  5  between  30  and  60  days  was  due 
to  reversed  surface  currents.  The  VAX  11/750  CPU  time  requirement  for  the 
WASP02  simulation  was  approximately  1  min. 

Comparison  with  WIFM-SAL 

39.  Two  different  injections  were  made  into  the  6:1  box  model  grid 
overlay  and  in  the  analogous  areas  of  the  WIFM-SAL  grid.  First,  a  25-mg/i 
spike  was  made  in  segment  683  of  the  WASP  model  (see  Figure  7)  and  the  six 
corresponding  cells  for  the  WIFM-SAL  simulation.  Second,  a  gradient-type  ini¬ 
tial  condition,  was  input  in  both  WASP  and  WIFM-SAL.  Initial  concentrations 
for  this  in  lection  are  tabulated  by  box  number  in  Table  4  and  shown  as  a  con¬ 
tour  plot  in  Figure  21.  For  a  direct  (1:1)  box  model  grid  overlay  of  the 
Mobile  Bay  area,  a  25-mg/i  dye  injection  was  made  at  WIFM-SAL  grid  lot. at  i  or. 

N  =  95,  M  =  8  (see  Figure  7). 

40.  Results  for  these  test  cases  are  presented  as  contour  plots  after 
4.5  days  of  simulation.  The  large  number  of  grid  points  and  box  segments  made 
plotting  individual  segment  concentration  histories  unwieldy.  Concent  ration 
contours  of  WIFM-SAL  results  after  4.5  days  of  simulation  for  the  25-mg/. 
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Figure  18.  DeGruv  lake  WASI'OJ  test  results  coinpa 
Ch-0L:AL-W2  results  after  30-day  simulation 
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Figure  19.  DeGray  Lake  WASP02  test  results  compared  with 
CE-QUAL-U2  results  alter  60-day  simulation 


Table  4 
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WASP  Gradient  Injection,  Initial  Concentrations,  Mississippi  Sound 


Cone 

Cone 

Cone 

Box  No. 

mg/ i 

Box  No. 

mg/  i 

Box  No. 

mg/ 

432 

0.1 

557 

5.0 

721 

5.0 

433 

0.1 

558 

5.0 

722 

1  .0 

434 

0.1 

559 

5.0 

723 

0.5 

435 

0.1 

580 

5,0 

724 

0.1 

436 

0.1 

581 

3.0 

747 

0.1 

437 

0.1 

582 

1.0 

748 

0.5 

438 

0.1 

583 

0.5 

749 

1  .0 

439 

0.1 

584 

0.1 

750 

3.0 

440 

0.1 

607 

0.1 

751 

3.0 

441 

0.1 

608 

0.5 

752 

3.0 

442 

0.1 

609 

1.0 

753 

3.0 

443 

0.1 

610 

3.0 

754 

3.0 

444 

0.1 

611 

5.0 

755 

3.0 

467 

0.1 

612 

7.0 

756 

3.0 

468 

0.5 

613 

7.0 

757 

1  .0 

469 

0.5 

614 

7.0 

758 

0.5 

470 

0.5 

615 

5.0 

7  59 

0.1 

471 

0.5 

616 

3.0 

782 

0.1 

472 

0.5 

617 

1  .0 

783 

0.5 

473 

0.5 

618 

0.5 

784 

i  .0 

474 

0.5 

619 

0.1 

785 

1.0 

475 

0.5 

642 

0.1 

786 

1  .0 

476 

0.5 

643 

0.5 

787 

1.0 

477 

0.5 

644 

1.0 

788 

1  .0 

478 

0.5 

645 

3.0 

789 

1 .0 

479 

0.  1 

646 

5.0 

790 

1.0 

502 

0.1 

647 

7.0 

791 

1  .0 

503 

0.5 

648 

9.0 

792 

1.0 

504 

1.0 

649 

7.0 

793 

0.5 

505 

1  .0 

650 

5.0 

794 

!  .0 

506 

1.0 

651 

5.0 

817 

0.1 

507 

1.0 

652 

1.0 

818 

0.5 

508 

1.0 

653 

0.5 

819 

0.5 

509 

1.0 

654 

0.1 

820 

0.5 

510 

1.0 

677 

0.1 

82! 

o .  :■ 

511 

1.0 

678 

0.5 

822 

0.5 
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1.0 

679 

1  .0 
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0.5 
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0.5 

680 

3.0 

824 

0.5 

514 

0.1 

681 

5.0 

825 

0.5 

537 

0.1 

682 

7.0 

8  26 

0.5 

538 

0.5 

683 

7.0 

82  7 

0.5 

539 

1.0 

684 

7.0 

828 

0.5 

540 

3.0 

685 

5.0 

829 

0.  1 

541 

3.0 

686 

3.0 

852 

0.1 

542 

3.0 

687 

1.0 

853 

0.1 

543 

3.0 
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0.5 

854 

0.1 

544 

3.0 
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0.1 

855 

0.  ! 

545 

3.0 
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0.1 
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Figure  21.  Mississippi  Sound,  WASP  gradient  injection 
initial  condition 


spike  injection  are  shown  in  Figure  22a.  WASP  results  for  tills  injection 
using  a  6-hr  time-step  are  shown  in  Figure  22b.  Clearly,  the  WASP  simulation 
is  overdispersive  in  comparison  with  the  WIFM-SAL  results.  Unlike  the 
CE-QUAL-W2  solution  scheme  (i.e.,  unlike  CE-QUAL-W2,  WIFM-SAI.  is  not  inher¬ 
ently  overdispersive),  the  flux-corrected  transport  scheme  in  WI FM-SAL  con¬ 
trols  numerical  dispersion.  The  spike  injection  into  the  Mobile  Bay  area, 
using  a  direct  box  model  overlay  on  the  WIFM-SAL  grid,  illustrates  that  for  a 
one-to-one  grid  overlay,  the  box  model  is  overdif fusive  in  comparison  with 
WIFM-SAL  transport.  Figure  23  demonstrates  this  using  time-steps  as  long  as  6 
and  12  hr  for  the  box  model  simulation.  Peak  concentrations  for  the  box  model 
simulation  are  twofold  to  threefold  lower  than  peak  values  maintained  by 
WIFM-SAL. 
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OCEPN  BOUNEiO^' 

a.  WIFM-SAL  results 

MOBILE  f - 1 


V.. 


DCfPN  RniJNOftRt 

b.  Box  model  results  for  1:6  overlay 

Figure  22.  Mississippi  Sound,  concentration  contours  after 
4.5-day  simulation  for  instantaneous  injection  in 

WASP  segment  683 
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a.  WIFM-SAL  results 


b.  Box  model  results  using  1:1  overlay  and 
6.0-hr  time-step 


c.  Box  model  results  using  1:1  overlay  and 
12.0-hr  time-step 


Figure  23.  Mobile  Bay  instantaneous 
injection  after  A. 5-day  simulation 


41.  For  most  water  quality  considerations,  a  gradient  oi  the  constitu¬ 
ents  rather  than  a  sharp  front  exists.  For  the  6:1  overlay  6-hr  time-step  and 
the  initial  conditions  described  in  Table  4,  results  are  shown  in  Figure  24. 
Although  the  overdispersive  character  of  the  box  model  is  still  discernible, 
the  box  model  more  closely  replicates  the  transport  for  the  gradient  injection 
condition  than  for  the  steep  front  condition. 
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PART  V:  SUMMARY  AMD  RECOMMENDATIONS 


42.  The  box  model  is  a  computationally  feasible  approach  for  long-term 
multidimensional  water  quality  modeling.  For  the  applications  described  in 
this  report,  the  box  model  required  under  15-min  CPU  time,  whereas  the  simula¬ 
tions  using  the  directly  linked  models,  CE-QUAL-W2  and  WIFM-SAL,  required  from 
3.0-  to  20.0-hr  VAX  11/750  CPU  time.  Computational  time  requirement  is  an 
important  consideration  for  a  complex  water  quality  application  requiring  a 
large  number  of  calibration  runs  and  long-term  (seasonal  and  annual) 
simulations. 

43.  The  box  model  is  numerically  diffusive.  An  increase  in  box  model 
segment  size  increases  numerical  diffusion.  Numerical  diffusion  decreases  as 
the  time-step  increases  until  numerical  instabilities  occur;  i.e.,  as  the 
Courant  number  approaches  one,  numerical  diffusion  approaches  zero. 

44.  Box  model  simulation  results  compare  well  with  results  obtained 
using  CE-QUAT.-W2.  Using  a  one-to-one  grid  overlay,  box  model  results  are 
nearly  identical  to  CE-QUAL-W2  results.  The  upwind  difference  scheme  in 
CE-QUAL-W2  is  also  a  diffusive  scheme  and  is  equivalent  to  the  box  nodel 
solution.  However,  even  using  a  one-to-one  grid  overlay  and  time-steps  up  to 
12  hr,  the  box  model  is  overdif fusive  in  contrast  to  WIFM-SAL,  which  has  a 
dif fusion-controlled  solution. 

45.  Neither  the  box  model  nor  CE-QUAL-W2  is  appropriate  for  simulating 
a  situation  where  movements  of  sharp  fronts  of  material,  such  as  spills,  must 
be  depicted  accurately.  Both  are,  at  present,  numerically  overdif fusive ;  i.e. 
numerical  diffusion  at  the  model  in  many  instances  exceeds  physical  diffusion 
in  the  system.  However,  for  situations  where  concentration  gradients  are  not 
steep,  these  models  can  adequately  represent  the  mass  transport  in  a  water 
body.  This  limitation  is  acceptable  for  many  water  quality  applications.  At 
the  time  of  this  publication,  both  models  were  being  modified  to  reduce 
numerical  diffusion. 

46.  Extensive  effort  to  accurately  determine  the  mass  dispersion  coef¬ 
ficients  for  a  multidimensional  modeling  study,  using  either  the  box  model  or 
CE-QUAL-W2  with  its  current  transport  algorithm,  appears  to  be  unnecessary, 
since  numerical  diffusion  in  these  models  is  typically  greater  than  the 
physical  dispersion  for  the  system.  However,  where  intertidal  averaging  is 


performed,  additional  dispersion  may  need  to  be  input,  and  calculation  of  this 
parameter  could  be  crucial. 

47.  Future  activities  should  focus  in  the  following  three  areas: 

a.  Linkage  of  the  box  model  with  other  hydrodynamic  models 
including  CELC3D,  a  three-dimensional  stretched  grid  model; 
CH2D,  the  two-dimensional  boundary-fitted  hydrodynamic  model; 
and  CH3D,  the  three-dimensional  boundary-fitted  hydrodynamic 
model . 

b.  Integration  of  a  dif fusion-controlled  solution  scheme  in  the 
box  model  and  CE~QUAL~W2  (see  paragraph  45), 

£.  Testing  of  the  use  of  intertidally  averaged  hydrodynamics  with 
the  box  model  and  development  of  methods  to  calculate 
advective  flows  and  dispersion  coefficients  for  the  box  model 
for  these  intertidal  cases  from  the  hydrodynamic  model  output. 
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